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INFLUENCE  OF  LEADING-EDGE  SUCTION  ON  LIFT-DRAG  RATIOS 
OF  WINGS  AT  SUPERSONIC  SPEEDS 
By  Clarence  B.  Cohen 


SUMMARY 

A  method  Based  on  linearized  theory  is  presented  for  calculat¬ 
ing  the  theoretical  suction  force  at  the  subsonic  leading  edges  of 
a  family  of  wings  at  supersonic  speeds.  The  method  is  used  to 
determine  the  optimum  sweepback  angles  for  the  tips  of  trapezoidal 
wings  and  to  determine  the  effect  of  curvature  of  the  tip  contour 
on  the  lift-drag  ratio  of  wing  regions  influenced  by  the  tip.  The 
effect  of  skin  friction  is  included.  The  possible  gain  in  lift- 
drag  ratio  from  proper  tip  design  of  trapezoidal  wings  increases 
as  the  sweepback  of  the  wing  is  increased.  Results  indicate  that 
appropriately  curved  tip  boundaries  will  give  higher  lift-drag 
ratios  in  the  region  affected  by  the  tip  than  the  best  trapezoidal 
wing. 


INTRODUCTION 

The  two-dimensional  theory  of  thin  airfoils  at  subsonic  speeds 
gives  a  resultant  force  normal  to  the  wing  surface  if  only  pressure 
forces  on  the  wing  surface  are  considered,  which  Implies  a  net  drag. 
For  this  theory,  however,  the  velocity  components  become  infinite 
at  the  leading  edge.  Evaluation  of  the  effect  of  this  singularity 
yields  a  suction  force  that  exactly  cancels  the  drag  and  thus  sat¬ 
isfies  the  momentum  requirement  that  the  resultant  force  vector  be 
normal  to  the  direction  of  flight  (reference  l). 

The  procedure  used  to  derive  the  suction  force  is  not 
restricted  to  two-dimensional  wings  and  indicates  that  such  a 
force  may  occur  along  any  edge  at  which  the  normal  component  of 
the  velocity  becomes  infinite.  In  the  linearized  theory  for  thin 
wings  flying  at  supersonic  speeds,  such  a  singularity  exists  along 
leading  edges  swept  behind  the  Mach  angle.  It  has  therefore  been 
proposed  by  Brown  (reference  2)  that  a  leading-edge  suction  force, 
which  tends  to  counteract  the  wave  drag,  exists  along  subsonic  v 
leading  edges  of  wings  flying  at  supersonic  speeds.  The  magnitude, 
of  this  force  is  calculated  in  reference  2  for  triangular  wings 
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swept  beyond  the  Mach  angle.  The  reduction  in  drag  predicted  by 
the  concept  of  leading-edge  suction  force  has  been  shown  by 
Vincenti  to  occur  experimentally  for  triangular  wings  swept  beyond 
the  Mach  angle. 

The  perturbation- velocity  components  required  to  evaluate  the 
local  leading-edge  suction  force  may  be  obtained  by  the  methods  of 
reference  3  for  the  class  of  wings  or  wing  tips  for  which  the  for¬ 
ward  Mach  line  from  the  subsonic  leading  edge  intersects  a  super¬ 
sonic  leading  edge.  An  analysis  of  the  suction  force  for  such 
wings  was  made  at  the  NACA  Cleveland  laboratory  in  1948  to  deter¬ 
mine  its  effect  on  the  lift-drag  ratio  obtainable  with  various  tip 
contours. 


SYMBOLS 

The  following  symbols  are  used  in  this  report: 
b  constant  defining  location  of  trailing  edge  of  wing 
^D,f  coefficient  of  drag  due  to  skin  friction 
C^  lift  coefficient 

C  pressure  coefficient 

Jr 

c  chord  of  wing 

D  drag  force 

F  suction  force  (in  flight  direction) 

G  quantity  defining  cpn 

k  constant  defining  angle  of  straight  wing  edge 

L  lift  force 

M  Mach  number  (free  stream) 

n  distance  normal  to  leading  edge 

q  dynamic  pressure,  ^pU2 
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E  distance  ratio,  u_u^  =  ^ 

S  integration  area 

s  distance  along  wing  edge 

U  free-stream  velocity 

u,  v  oblique  coordinates 

V  distance  (in  v-direction)  between  subsonic  and  supersonic 
leading  edges 

x,  y  Cartesian  coordinates 
a  angle  of  attack,  radians 

3  cotangent  of  Mach  single,  vM2-l 

e  semiapex  angle  of  triangular  wing 

0  angle  between  wing  edge  and  flight  direction  (positive 
counterclockwise),  . ^ 

li  Mach  angle,  sin"1  ^ 

p  mass  density  of  air 

cp  perturbation-velocity  potential 

Subscripts: 

1  supersonic  leading  edge 

2  subsonic  leading  edge 

f  friction 

n  indicates  direction  normal  to  wing  leading  edge 
t  tip  point  of  wing  (at  greatest  span) 

w  wave 

x,  y  indicates  velocity  components  in  x,  y  directions 
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x' 


DERIVATION  OF  EXPRESSION  FOR  LEADING-EDGE  SUCTION  FORCE 

In  subsonic  airfoil  theory,  the  form  of  the  equation  of  the 
velocity  component  normal  to  the  leading  edge  (reference  l)  is 


*n 

n->0 


G 


(i) 


where  cpn  =  ^  is  the  perturbation  velocity  in  the  n  direction 

(fig.  l).  Thus  cp^  approaches  infinity  as  n  approaches  zero. 

Analysis  of  this  singularity  results  in  a  suction  force  per  unit 
length,  which  can  be  expressed  in  differential  form  as  (refer¬ 
ence  1) 


dF  „ 

=  *pG  (2) 

where  s  is  distance  along  the  wing  edge  (fig.  1) . 

Equations  (l)  and  (2)  may  be  rewritten,  as  in  reference  2,  to 
include  the  effects  of  compressibility. 


and 


(la) 


(2a) 


where  is  the  component  of  the  Mach  number  normal  to  the  edge. 

For  subsonic  two-dimensional  wings,  the  factor  G  is  such 
that  the  force  defined  by  equation  (2a)  is  just  sufficient  to 
cancel  the  drag  resulting  from  the  flat-plate  analysis.  For 
supersonic  wings,  the  factor  G  must  be  determined  from  the 
appropriate  value  of  cpn  at  the  subsonic  leading  edge.  From 

reference  3,  the  perturbation- velocity  components  in  the  x  and  y 
directions  near  the  leading  edge  can  be  written  (appendix  A)  as 
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The  coordinate  system  used  herein  is  shown  in  figure  1.  The 
relations  between  the  oblique  and  Cartesian  coordinates  are  given 
by 


Equations  (3)  may  be  combined  to  obtain  the  normal  velocity 
component  near  the  edge 


u-u2  u-u2 


The  geometrical  relation  between  E  and  n  as  n->0  is 
shown  in  figure  2,  where 

n  =  PQ  sin  (0-Ha)  =  ^  sin  (0-Ha)  (6) 

The  following  expressions  are  derived  from  the  relations 
between  the  oblique  and  Cartesian  coordinates: 
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sin  6  = - - - 

VT 


cos  9  =  — - — — 

•Jk- 

With  the  aid  of  equations  (6)  and  (7),  equations  (la)  and  (5) 
may  he  solved  for  G2. 


p 

This  value  of  G  substituted  into  equation  (2a)  yields  the 
normal  force  per  unit  length  along  the  wing  edge. 


^n  _  4PU2a2V 
ds  itM  N  dv 


With  the  use  of  the  relation 


dF  dF  ds  ds 

dv  ds  dv  ds  dv  8 


(9) 


u 


TOO! 
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there  results 


ar  _  ioi^g^v 

dT' 

The  leading-edge  suction  force  can  now  be  expressed,  in  inte¬ 
gral  form  as 

'■‘ff-  /•(> 


(  du2\  m* 

The  quantity  ll  -  in  equation  (11)  is  given  in 

figure  3  as  a  function  of  0  for  three  values  of  3.  With  the 
aid  of  curves  such  as  these  and  a  drawing  of  the  wing  in  question, 
equation  (ll)  may  he  evaluated  by  a  numerical  summation. 

The  force  defined  by  equation  (ll)  is  obtained  when  the  radius 
of  curvature  of  the  leading  edge  is  vanishingly  small  and  when  the 
velocity  components  normal  to  the  edge  become  infinite  at  the  edge. 
In  practice,  however,  the  infinite  velocity  components  predicted  by 
linearized  theory  correspond  to  a  high  negative  pressure  coefficient 
near  the  edge.  By  rounding  the  subsonic  leading  edge  to  increase 
the  projected  area  on  which  the  lower  pressure  can  act,  it  is 
possible  to  reduce  the  total  drag  of  the  wing.  As  the  edge  radius 
is  increased,  an  optimum  curvature  that  corresponds  to  a  maximiza¬ 
tion  of  the  product  of  area  times  suction  force  per  unit  area, 
should  be  obtained.  It  has  been  shown  experimentally  that  round¬ 
ing  the  leading  edge  of  triangular  wings  swept  beyond  the  Mach 
angle  increases  the  lift-drag  ratio,  but  that  continued  rounding 
beyond  a  certain  value  produces  no  further  gains. 


APPLICATION  OF  LEADING-EDGE  SUCTION-FOKCE  EQUATION 

For  a  given  main  surface  (fig.  1),  the  optimum  wing  is  one 
whose  tip  has  the  greatest  lift-drag  ratio  L/D.  If  the  tip  has 

1 

no  subsonic  edge,  this  ratio  is  simply  i  =  -  -  -  =  — - — 

D  Dw  +  Df  .  £f 

J-  »  TV 


where 
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D  is  the  wave  drag  and  the  skin-friction  drag.  If  the  tip 

W  I 

has  a  subsonic  leading  edge  its  lift-drag  ratio  becomes 


L  _  L 

D  D  +  D-  -  F 

W  X 


(12) 


where  F  is  the  suction  force  at  the  subsonic  leading  edge.  The 

/Df  \ 

tip  with  the  minimum  value  of  I  g —  •g-J  therefore  has  the  highest 

lift-drag  ratio.  If  frictionless  flow  is  as  sinned,  the  param¬ 
eter  f/Dw  should  be  maximized. 

As  an  example  of  the  procedure  for  obtaining  wings  of  optimum 

(%  F  \ 

plan  boundaries,  the  quantity  l  -  ■g-  1  will  be  derived  for  a 

wing  with  a  straight-line  supersonic  leading  edge  v^=  k^u.  The 
quantity  V  in  equation  (ll)  then  becomes  v  -  k^Ug  (fig.  4)  and 
the  suction  force  is 

f  ■  ^  f  *  <',-ki"2)  (x  •  ?)  \S- dT  (i3) 

Jo 

where  v^  is  the  value  of  v  at  the  extreme  tip  and  u  =  u^v) 
is  the  equation  of  the  tip  edge. 


The  wave  drag  for  a  wing  tip  of  this  type  is  (appendix  B) 


4pP2»2 

V.  __w2 


tMk- 


J~h. 


1  +  ^  ^-2  a^(v- k^)  [b-Cv-k^Ug)] 


/b-(v- k-jUg) 

+  (kr1)b  tan  V  (v-k^)- 


(14) 


and  the  friction  drag  is 
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^  "  *  (l~kz)  klA/'klk2 
Dw  (3kx  -  2k1k2  -  l) 

where 

-  00  <  k-^  <  0 

0  <  kg  <  1 

and  and  kg  determine  the  sweep  angles  of  the  leading  and 

trailing  edges  0^  respectively.  The  angles  corre¬ 

sponding  to  several  values  of  the  constants  are: 

ki  ei 
0  -n 
-1  «/2 
-00  |J. 

The  equations  relating  the  angles  and  the  constants  axe 


(18) 


As  an  indication  of  the  variation  of  leading-edge  suction 
force  for  this  type  of  wing,  the  effect  of  02  upon  F/D^,  as 

calculated  from  equations  (17)  and  (18)  with  CD  f  =  0,  is  shown 
in  figure  5  for  0^  =  65°  (k^  =  -2.75)  at  M  =  ,^2.  An  optimum 

value  of  0g  occurs  at  about  20°. 

In  order  to  obtain  an  expression  for  optimum  values  of  02 
(or  k2)  as  a  function  of  k^,  equation  (17)  must  be  differ¬ 
entiated  and  set  equal  to  zero.  When  differentiated  with  respect 
to  k2,  the  equality  yields 


4 
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2kg2  kx  +  k2  (3-7^)  +  3kx  -  1  +  JtPC^-f  /k^  =  0  (19) 

or  y 

Differentiation  with  respect  to  k^  yields  no  true  nwMriimnn  for  the 
type  of  wing  treated.  The  variation  of  the  optimum  value  of  k2 
with  k-, ,  as  expressed  by  equation  (19),  is  shown  in  figure  6  for 

PCDf 

several  values  of  — x2— . 

or 

A  wing  with  ©^  =  65°  and  three  alternate  tip  curves  A,  B, 

and  C  passing  through  points  0  and  t  is  shown  in  figure  7.  Tip 
curve  B  is  swept  so  that  ©g  =  20°  (the  optimum  for  this  trape¬ 
zoidal  wing  from  fig.  5).  Curves  A  and  C  are  arbitrary  curves, 
symmetrical  about  B,  drawn  through  0  and  t. 

The  following  values  of  f/b v  were  obtained  by  numerical 

integration  of  equations  (13)  and  (14)  for  these  tips;  the  result¬ 
ing  effect  on  lift-drag  parameter  aL/D  for  frictionless  flow  is 
also  shown: 


Wing 

»/>v 

oL/D 

(frictionless 

flow) 

A 

0.39 

1.64 

B 

.31 

1.45 

C 

.20 

1.25 

These  results  indicate  that  the  lift-drag  ratio  may  be  Increased  by 
the  use  of  appropriately  curved  wing  boundaries  at  supersonic 
speeds . 

As  an  indication  of  the  magnitudes  of  F/Dy  theoretically 
obtainable  with  other  type  plan  forms,  the  variation  of  F/hy 

with  apex  angle  for  a  triangular  wing  is  shown  in  figure  8.  This 
relation  is  based  on  the  equations  in  reference  2. 

The  variation  of  maximum  lift-drag  ratio  with  k-j^  is  shown 

in  figure  9  for  trapezoidal  tips.  These  lift-drag  ratios  were 
calculated  from  the  optimum  values  of  k2  shown  in  figure  6. 
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The  lift- drag  ratio  for  the  main  surface  of  a  trapezoidal  wing  is 
also  shown  for  comparison.  These  curves  show  that  the  maximum 
lift-drag  ratio  of  sweptback  trapezoidal  wings  is  obtained,  not 
when  the  wing  is  cut  off  at  the  inner  Mach  line  (main- surface 
curves),  but  rather  when  the  optimum  tip  angle,  based  on  considera¬ 
tion  of  leading-edge  suction,  is  chosen.  (This  result  held  for  all 


Pcn  f 

values  of  the  friction  parameter  — shown. )  For  sweptforward 


a 


trapezoidal  wings,  however,  figure  9  indicates  that  the  optimum 
lift-drag  ratio  L/D  is  obtained  when  the  tip  is  cut  back  at  the 
inner  Mach  angle. 


The  lift-drag  ratio  is  shown  in  figure  10  as  a  function  of 
angle  of  attack  for  the  wing  tips  previously  mentioned,  for  the 
main  surface,  and  for  a  tip  cut  off  in  the  flow  direction  (tip  D). 
The  calculations  were  made  for  M  =  */2  and  C-^f  =  0.005.  This 

figure  indicates  that  appropriately  curved  plan  forms  give  more 
favorable  lift-drag  ratios  than  the  optimum  trapezoidal  wing.  For 
wings  of  low  aspect  ratio,  the  effect  of  the  tip  design  becomes 
relatively  more  important;  a  complete  study  must  therefore  include 
the  effect  of  aspect  ratio.  From  the  standpoint  of  high  lift-drag 
ratio,  a  tip  such  as  D  (fig.  10)  is  poor  for  any  aspect  ratio 
inasmuch  as  a  higher  lift-drag  ratio  can  be  obtained  merely  by 
cutting  off  the  wing  at  the  Mach  line. 


CONCLUSIONS 

The  following  conclusions  result  from  a  method  for  evaluation 
of  the  theoretical  leading-edge  suction  force  on  a  type  of  wing 
for  which  the  forward  Mach  lines  from  the  subsonic  leading  edge 
intersect  a  supersonic  leading  edge: 

1.  The  optimum  lift- drag  ratio  of  sweptback  trapezoidal  wings 
is  obtained  when  the  optimum  tip  angle,  based  on  consideration  of 
leading-edge  suction,  is  chosen.  For  sweptforward  trapezoidal 
wings,  however,  the  optimum  lift-drag  ratio  is  obtained  when  the 
tip  is  cut  back  at  the  Mach  angle.  The  possible  gain  in  lift-drag 
ratio  from  proper  tip  design  of  trapezoidal  wings  increases  as  the 
sweepback  of  the  wing  is  increased. 

2.  Higher  lift-drag  ratios  may  be  obtained  with  wing  tips 
whose  edges  are  appropriately  curved  than  may  be  obtained  with 
straight-line  wing  tips. 


Lewis  Flight  Propulsion  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Cleveland,  Ohio,  July  15,  1948. 
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APPENDIX  A 


PERTURBATION- VELOCITY  COMPONENTS 
NEAP  SUBSONIC  LEADING  EDGE 

Sep 

The  x- component  of  the  perturbation  velocity  cp  =  —  where 

x  ox 

cp  is  the  perturbation-velocity  potential,  is  given  in  reference  3 
as 


Sep  _  Uq 
Sx 


dv'-du' 


(u— U *  )  (  V— V  *  ) 


^(u-u' )(v-v*) 


(Al) 


where  the  primes  denote  the  variables  of  integration. 


The  integrals  are  evaluated  along  the  supersonic  leading  edge 
from  A  to  B  and  along  the  Mach  line  from  B  to  Q.  When  point  P 
approaches  the  subsonic  leading  edge  (5q  (n-*0),  points  A  and  B 

converge  and  equation  (Al)  reduces  to 
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n-X> 


_ dv1 _ 

V(u-U* ) (v-v* ) 


(A2) 


which  upon  integration 
expression  for  cp^.  is 


becomes  the  first  of  equations 
obtained  similarly. 


(3). 


The 
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APPENDIX  B 


DRAG  FORCE 

In  accordance  with  linearized  theory,  the  Wave  drag  of  a  thin 
flat  plate  is 

D^  =  aL  =  aCLqS  (Bl) 

where  the  lift  coefficient  can  he  expressed  as 


For  a  wing  of  the  type  shown  in  figure  4,  reference  3  shows 
that  the  pressure  coefficient  can  he  expressed  as 


C  =  - 


fl  -  ^  /v~klu2  +  1~kl 


[k-^ug-u) 

v"klu2 


A  combination  of  equations  (Bl)  to  (B3)  yields 


which  upon  integration  with  respect  to  u  becomes  equation  (14). 
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APPENDIX  C 


SUCTION  FORCE  AND  DRAG  FOR  A  TRAPEZOIDAL  WING 
The  equation  of  the  tip  of  a  trapezoidal  wing  is 

u^v)  =  kgv  (Cl) 


By  the  use  of  this  equation,  if  the  value  of  is  known,  the 

wave  drag,  friction  drag,  and  suction  force  can  be  obtained  from 
equations  (13)  to  (15).  For  the  trapezoidal  wing 

-  ,  L_c_\  *  x>2  -  4kl  (02) 

t  \l-^2)  2p 

The  use  of  equations  (Cl)  and  (C2)  in  equations  (13)  to  (15)  yields 


*  -  py  ph  * *>»  -  m  <«> 


D  puVc2  j  [m2 (^l  +  l)2  -  4kJ  (3k!  -  2kxk2  -  1)1 

8P2  |  (l-kikg)  klN/^[ 


(C4) 


pU2c2CD>f  [m2^  +  l)2  -  4k  J 

Df  =  eT-2  k^vD 


(C5) 


By  subtracting  equation  (C3)  from  (C5)  and  by  dividing  the  result 
by ‘equation  (C4),  equation  (17)  is  obtained. 
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Wing  boundary 


v  -  vx(u) 


R  a-5£-  JL. 
PQ  u-u2 


4  -1  1 

sin  A 


n  a  X  sin  (0+M-) 


(e+n] 


P(u,v) 


Figure  2.  -  Geometric  relation  between  R  and  n. 


Figure  4.  -  Wing  with  straight  supersonic  leading  edge  and 
generalized  subsonic  leading  edge  up(v). 
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Figure  6.  -  Variation  of  kg  with  for  maximum  lift-drag  ratio  of  trapezoidal 

wing  tips. 
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Figure  9*  -  Maximum  lift-drag  ratio  for  regions  of  trapezoidal  wings  when  effects 
of  leading-edge  suction  and  friction  drag  are  considered. 
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d/1  *o 


Angle  of  attack,  a,  deg 

Figure  10.  -  Variation  of  lift-drag  ratio  with  angle  of  attack  at  M  =  42*  and 

CL\  a  *  0*005. 


